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This paper contains a perturbation theory which is applicable to the 
scattering losses suffered by guided modes of a dielectric slab waveguide as a 
consequence of imperfections of the waveguide wall. The development of 
the theory occupies the bulk of the paper. Numerical results appear in 
Sections VI and VIII to which a reader less interested in the theory is 
referred. 

The theory allows us to conclude that random deviations of the waveguide 
wall in the order of 1 percent, for guides designed to guide an optical wave 
of X = In wavelength, can cause scattering losses of 10 percent per centi- 
meter or O.46 dB per centimeter. A systematic sinusoidal deviation of the 
waveguide wall can cause total exchange of energy from the lowest order to 
the first order guided mode in a distance of approximately 1 cm if the ampli- 
tude of the sinusoidal deviation from perfect straightness is only 0.5 percent 
of the thickness of the guide. An rms deviation of one of the waveguide walls 
of 9 A causes a radiation loss of 10 dB per kilometer (index difference 
1 percent, guide width 2.5 n). 

I. INTRODUCTION 

The problem of how to transmit laser light over large distances or 
carry it short distances inside the laboratory has renewed the interest 
in dielectric waveguides. 1-5 Such waveguides usually used in the form 
of clad fibers or as strips of a medium of larger dielectric constant 
embedded in another dielectric medium are capable, in principle, of 
guiding electromagnetic radiation. By proper dimensioning, a dielectric 
waveguide can be made to transmit only one guided mode. In this 
respect mode guidance by dielectric waveguides resembles mode guid- 
ance by hollow metallic waveguides. Hollow metallic tubes can be 
constructed to allow only one mode to propagate so that mode conver- 
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sion (except for conversion to the reflected dominant mode) becomes 
impossible. Such truly single mode operation is impossible for dielectric 
waveguides since these guides can always lose electromagnetic energy 
to the continuous spectrum of unguided modes. 

The possible solutions of Maxwell's equations for a dielectric wave- 
guide consist of a discrete spectrum of a finite number of guided modes 
plus a continuum of waveguide modes. 6 The guided modes have field 
configurations which concentrate the electromagnetic energy inside 
and in the immediate vicinity of the structure. The continuum of un- 
guided modes extends to infinite distances from the waveguide and 
consists of a superposition of incident and reflected waves. A convenient 
way of visualizing the physical significance of the continuum of unguided 
modes is as follows. If a plane wave is incident on the dielectric wave- 
guide at an arbitrary angle, part of it penetrates the dielectric structure 
while some portion is reflected. The resulting superposition field of 
incident and reflected waves satisfies Maxwell's equations and the 
boundary conditions at the dielectric waveguide and as such can be 
viewed as a mode of the structure, but the energy of this mode is not 
concentrated near the waveguide and there are no specific restrictions 
on the projection of the propagation vector in the direction of the guide 
axis. 

A perfect dielectric waveguide can transmit any of its guided modes 
without converting energy to any of the other possible guided modes 
or to the continuous spectrum. But any imperfection of the guide, such 
as a local change of its index of refraction or a deviation from perfect 
straightness or an imperfection of the interface between two regions 
with different index of refraction, couples the particular guided mode to 
all other guided modes as well as to all the modes of the unguided con- 
tinuum. Imperfections of this type are unavoidable. They transfer 
energy from the desired guided mode to unwanted guided modes and 
the radiation field of the continuum of unguided modes, thus increasing 
the loss of the desired guided mode. 

This paper gives a simple, approximate theory of the losses of di- 
electric waveguides, caused by imperfections of the boundary between 
the inner region of higher dielectric index and the surrounding outer 
region of the dielectric waveguide. Even though the method of analysis 
used here can be used to describe any arbitrary dielectric waveguide, 
we limit the discussion to a simple case. We describe the effects of mode 
conversion for a dielectric slab surrounded by vacuum, assuming for 
simplicity, that there is no variation of the dimensions or properties of 
the rod as well as the field distribution in one co-ordinate direction. The 
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restriction of demanding d/d v = for one of the co-ordinates y is no 
limitation on the method of analysis but is imposed strictly for con- 
venience. It simplifies the analysis considerably without drastically 
changing the conclusions. The tolerance requirements based on our 
analysis are rather stringent. They show the order of magnitude of the 
losses which can be expected from deviations from perfect geometry. 
Additional variations in the direction considered perfect in this paper 
is unlikely to improve any of the loss predictions. 

II. TE MODES OF A DIELECTRIC SLAB 

Let us consider the transverse electric modes of the dielectric slab 
of Fig. 1. True to the simplifying assumption discussed in Section I, 
we assume 






(1) 



with y being the co-ordinate perpendicular to the x and z directions, but 
parallel to the slab. The only nonvanishing field components are E v , 
H t , and H. . 

Leaving the z and time dependence 

«*<—'•> (2) 

understood, we obtain the following modes of the ideal structure as a 
solution of Maxwell's equations satisfying the boundary conditions. 
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Fig. 1 — Geometry of a dielectric slab waveguide. 
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2.1 Even Guided Modes 
For even guided modes 

8„ = A <r) cos kx for | x | g d, (3a) 

S„ = .4 ( " cosKc/e- 7 "-" for .r ^ d, (3b) 

CO/J OZ 

K. = -ff«, (5) 

co^i o.r 

The field component 8„ satisfies the wave equation 
<^ + ^ + ni;/c 0. 

d.E dz 

The value of the index of refraction n is different inside and outside of 
the dielectric slab. For simplicity, we assume 

n = 1 for | x | > d. (7) 

The other constants are related as follows 

k 2 = <o 2 e /xo , (8) 

k = (n 2 k 2 - a )», (9) 

7 = (/3 2 " fc 9 ) 1 - (10) 

The propagation constant # is obtained as a solution of the eigenvalue 
equation 

tan Kd = -• (11) 

K 

The mode amplitude A can be expressed in terms of the power P carried 
by the mode. 

P = \ Re f (- 8„0C X *) dx = — f | 8„ | 2 tf.r. (12) 

P is the power per unit length (unit length in y-direction) flowing along 
the 2-axis. We obtain for the amplitude coefficient 

r = -^?. as) 

#* + * 
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2.2 Even Modes of the Continuum 

The continuum of unguided modes of even sj^mmetry is given by the 
equations: 

8„ = B ie) cos trx for \x \ ^ d, (14a) 

8„ = C {t) e<" + D M e~ ip ' for x ^ d. (14b) 

The other field components follow again from equations (4) and (5) 
and £„ is a solution of equation (6). The constants are related to each 
other by the equations 

o- = (n 2 k°- - t?)\ (15) 

p = (k 2 - /3 2 )>. (16) 

The radial propagation constant p can assume all values from to <» . 
The continuous mode spectrum starts at = k and continuous to /3 = 
at which point we have p = k. Larger values of p are obtained for 
imaginary values of /3 corresponding to modes of the continuum exhibit- 
ing a cutoff behavior. 

The boundary conditions do not lead to an eigenvalue equation for 
/3 but they determine C (e) and D (c) in relation to B U) . 

C le) = %B (e) e- ipd (cos ad + i- sin ad) , (17) 

D (e) = C {c) *, (18) 

(the asterisk indicates the complex conjugate quantity). 

The normalization of the modes of the continuum involves a 6-func- 
tion. Instead of equation (12) we use 

P o(p - p') = -2- [ S,( P )S*(p') dx. (19) 

cop Jo 

With this normalization we get 



7r/3( cos 2 ad H — 2 sin 2 ad 



2 .3 Odd Guided Modes 



In a maimer similar to that for obtaining the preceding equations we 
obtain the equations for the odd guided modes 

S y = A m 8inicx for x ^ d, (21a) 
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8, = A (0) am K de" ,l '-* ) for x ^ d. (21b) 

Equations (4) through (10) apply to the odd modes unaltered. The 
eigenvalue equation is given by 

(22) 



P. (23) 



2 .4 Odd Modes of the Continuum 

As in Section 2.3 we obtain the equations for the odd modes of the 
continuum 



tan nd 


K 


and the mode normalization is 




A (0) ' = 


2wn 

7 



8„ = B (0) sin ax for | x | ^ d, 


(24a) 


S„ = C w e ipx + D (0) e-' pi for x ^ d, 


(24b) 


C (o> = i£<°> e -->4 sin a d - i 1 cos (rdj , 


(25) 


7) (0) = r*^ 5 * 


(26) 


^(o,. 2co M P 


(27) 


~ ( a* \ 


ir/3(sin 2 trd H 2 COS 2 trrf ) 





All these modes are orthogonal to one another. The even modes are 
orthogonal to all the odd modes, the guided modes are orthogonal to 
all the modes of the continuum, and all guided modes as well as all 
modes of the continuum are orthogonal among each other. The or- 
thogonality of the modes of the continuum among each other was 
already expressed by equation (19). Labeling the discrete modes by 
indices and dropping the vector component label y we can express the 
orthogonality of the discrete modes by the equation 

P S nm =^/°° m S » 8 *^- (28) 

III. MODE COUPLING CAUSED BY IMPERFECTIONS 

We want to study the losses which the lowest order guided mode 
suffers because of imperfections of the waveguide wall. A dielectric 
waveguide with wall imperfections is shown in Fig. 2. 
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Fig. 2 — Dielectric slab waveguide with wall distortions. 

The waveguide with wall imperfections is mathematically described 
by a refractive index distribution 



n 2 (x, z) = n 2 (x, z) + An 2 {x, z). 
The index distribution 

\x\ < d 



n (x,z) = 



x I > d 



(29) 



(30) 



describes the ideal dielectric waveguide whose TE modes were given in 
the Section II. The additional term An 2 describes how the guide deviates 
from its perfect shape. Consider a deviation shown in Fig. 3. The 
corresponding distribution An 2 is (n„ = index of refraction of the di- 
electric material of the guide) 

\x < d if d < j(z) 

\x < Kz) if d > m 

n] - 1 d < x < f(z) if d < j(z) 

-(n] - 1) /(z) < x < d if d > j(z) 
'x > j(z) if d < f(z) 



An = 



(3D 



0- 



[x > d if d > f(z) 
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I \x = f(z) 



Fig. 3 — Illustration of the wall distortion function j(z). 



The field distribution E v of this waveguide is a solution of 
' r/ ' 1 " + ^§* + (nl + An 2 )k 2 E u = 



dx* 



dz 



(32) 



with H x and H t given by equations (4) and (5) . The modes of the perfect 
waveguide form a complete orthogonal set for all TE modes with no 
variation in the ^-direction. It is, therefore, possible to express any 
field distribution on the waveguide with imperfect walls by the expansion 



E, 



= E C„(*)s n + S f 9(p, z)e(p) dp. 

n JO 



(33) 



The first summation extends over all even and odd modes of the discrete 
spectrum of guided modes. The integral extends over all modes of the 
continuum, and the summation sign in front of the integral indicates 
summation over even and odd modes. The expansion coefficients C n 
and g(p) are unknown functions of z. 

To obtain a coupled system of differential equations for the expansion 
coefficients we substitute equation (33) into equation (32). Multiplying 
the resulting equation by 

2uu bm ' 

integrating over x from -co to + °° , and using the orthogonality rela- 
tions and the fact that £„ and &(p) are the (discrete and continuous) 
modes of the perfect guide leads to 



d 2 C m . a dC m „ . N 



(34) 



i." * 
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with 

FM = ~l^p [e CM /_" S* An 2 S„ dx 

+ E /f (/ pf/(p> 2 ) /_" s * Aw * 6 <p> dx J • 



(35) 



Similarly multiplying by 



if" S*(p') 



leads to 



with 



«^ - 2#' -^ = G(p'») (36) 



r;(p '- s) = ~£^ [? c " w .C 8 * (p0 An2£ - dx 

+ E r dpg( P ,z) f g*(p') An 2 8(p)dx • (37) 

No n-label on the power term P is necessary since we assume that all 
the normal modes are normalized to the same amount of power. The 
actual power carried by each mode relative to the power of the other 
modes is given by the C„ coefficients. Solutions of equations (34) and (35) 
with appropriate initial conditions provide us with exact solutions of 
the imperfect waveguide. It is interesting to note that this method of 
solution does not require the consideration of boundary conditions. 
The normal modes S„ and £(p) were assumed to have the time and 
z-dependence of equation (2); this means they represent waves traveling 
in the positive z-direction. However, the solutions of equations (34) and 
(36) introduce waves traveling in positive as well as negative z-direction. 
To see this, let us assume that An 2 = so that F„(z) = 0. The equation 

d °' C " 2i(i n — n = (38) 



dz ™ dz 

has the solution 

C n (z) = A + Be 2i0 »' (39) 

with constant A and B. The product of A with S„ results in a wave 
traveling in the positive z-direction but the product of B exp (2i'/3„z) 
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with 8 n results in a wave traveling in the negative z-direction. So, even 
though we started out with waves traveling in the positive z-direction 
the expansion (33) contains partial waves traveling in positive as well 
as negative z-direction. 

For the purpose of obtaining perturbation solutions of equations 
(34) and (36), an integral form of these equations is more useful. Treat- 
ing equations (34) and (36) as inhomogeneous differential equations, 
we can immediately write the following integral equations 

C m = A m + BJ^' + =%-(' [e 2iPmil - n - l]F.(f) df , (40) 

g( P ', z) = C( P ') + JXaV"' - + ^ [ V' 1 "" ~ W^'» f) d * ■ 

(41) 

It is important to know which part of equations (40) and (41) is as- 
sociated with waves traveling in the positive or negative z-direction. 
Therefore, we introduce the notation. 

C m = Ci +) + Ci _) (42) 



with 



C 



(+> 



» = A ~ - W m /.' FM) dr ' (43) 

Cl-\z) = {j5 m + -L- £ e- UM Fjf) dry*-'. (44) 

The superscript (+) indicates the coefficient which after substitution 
into equation (33) produces waves traveling in positive z-direction, 
while (— ) indicates the part which produces waves traveling in negative 
z-direction. A similar notation and resulting equations is used for g(p', z); 
however, the corresponding equations are obvious and are therefore 
omitted. 

The constants A m , B m , and so on, occurrng in equations (43), (44), 
and the corresponding equations for g(p', z) must be determined from 
initial conditions. We always assume that the lowest order guided mode 
is incident on the imperfect waveguide at z = 0. Using the subscript 
for this incident mode we get immediately from equation (43) 

Cl +) =0 for m^O at 2 = 
or 

A m = for m ^ 0, (45) 
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but 

A = 1. (46) 

We imagine that at z — L the waveguide is connected to a perfect guide 
so that at that point there are no waves traveling in negative z-direction. 
This leads to the condition 

B ~ " -5S- I " e- 2iM F m (t) dt (47) 

for all values of m. The power loss AP of the incident mode due to mode 
conversion is given by 

¥ - E [| c^\L) | 2 + i cr'(o) I 2 ] 

* n-1 

+ Zf[| <? (+> (p, L) | 2 + | ^"'(p, 0) | 2 ] dp. (48) 

•'o 

Equation (48) states that the total power lost by mode conversion from 
the incident mode escapes at z = L in spurious modes traveling in posi- 
tion z-direction and at z = in spurious modes traveling in negative 
z-direction. The factor P is the normalized power factor of equations 
(12) and (19); it is the power incident in mode 0. Notice that because 
of equations (45) and (47) only the integral terms of equations (43) and 
(44) (taken from z = to z = L) enter into equation (48). 

The integral equations (43) and (44) can only be solved approxi- 
mately. We perform first order perturbation theory by using C m (0) 
instead of C m (z) and g{p, 0) instead of g(p, z) in equations (35) and (37). 
Furthermore, we realize that CJ, -> (0) for all m is a quantity of first 
order and will therefore be neglected in equations (35) and (37). The 
same is true for C l m +) (0) with m 9* 0. In the spirit of first order perturba- 
tion theory we use therefore 

C m = 8 0m (49) 

and 

g(p) = (50) 

in equations (35) and (37). 

The perturbation theory is feasible not only when n] — 1 <3C 1 but 
also when n 2 — 1 is arbitrarily large but the geometrical deviation of 
the guide walls from perfect straightness is slight. In either case we 
obtain from equations (35) and (37) the simple approximations 
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F,„(z) = -§^>« ' l)IEfW - d\8 (d,z)&*(d,*) 

- [h(z) + d\&o(-d,z)&*(-d,z)), (51) 
G(p,z) = -£^p(ni " Dm) ~ m*(P,d,z)Go(d } z) 

- [h(z) + d] S*( Pj -(/, 2) S„(-r/, z)}. (52) 

The function /(z) describes the dielectric-air interface in the vicinity 
of x = d, while h(z) describes it near z = —d. We assumed that j(z) 
and h(z) depart so little from x = d and x = —d that the functions 
S(x, z) could be replaced by &(±d, z). 

IV. EVALUATION OF THE SPURIOUS MODE AMPLITUDES 

We begin the discussion of the consequences of our scattering theory 
by calculating the coefficients Cl +) and <7 < + ) . We obtain [from equations 
(43) and (51) with the help of equations (3a) and (13) for the even 
modes] the following 

CIV(L) = g (n! - 1) cos^cos^ ^ _ ^ (g3) 



2iV " "[(w+;ia^+5 j 

The coefficients tp m and \f/ m are defined by 

*>„ = 7 C U(z) ~ d]e- iifi -' mU dz (54) 

L/ .'o 

and 

*. - 7 r Wfl + flr ,cfc -* ,>- <&■ (55) 

These are the Fourier coefficients of the functions /(z) - d and /i(z) + d 
which are expanded in a domain 

S z £ L. 

The amplitude of the with even mode depends on the Fourier components 
of the wall function whose "spatial frequency" T is 

r„, = 2t = Po - p. . (56) 

The corresponding expression for the even modes of the continuous 
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spectrum is: 

(+)/ ,, Ur , 2 cosk„(/ cos aciyip) - MQ] 

^/ + ^(cos 2 o-r/ + ^sin 2 a ( / 

with [/3 = /3(p) see equation (1G)] 

*G8) = \ [ U(A - d\e- il0 '- fiu dz, 

W) = l f Q [*0 + d]e- iif, '- pu dz. 
The corresponding expressions for the odd modes are 



(57) 



(58) 
(59) 



Lk 



CffCL) = T5" (»' - 1) r 



cos Kod sin K„d 



2?: 



7o/\ 7n/J 



j fa. + *„), (60) 



( + ) . . Lk 2 , 2 cos K d sin o-rifcffi) + ^Q3)] 

00 (p, Ir) = ^TTT (n, - 1) f7 TT-7 2 

( ^ (/Sod + ^Wsin 2 <rtf + 2, cos 2 <rd 

_\ To/ \ P 

(61) 

The Fourier coefficients <p and ^ are given by equations (54), (55), 
(58), and (59) except that /3„ and /3 are now the propagation constants 
of the odd modes. 

The corresponding expressions for C <_) and g { ~ } are obtained by 
replacing /3 m with —0 m and with — in equations (54), (55), (56), (58), 
and (59)'. 

V. SINUSOIDAL WALL DEFLECTIONS 

As a specific example, let us assume that the wall imperfections have 
sinusoidal shape. Then 



and 



j(z) — d — a sin dz 



h(z) + d = —a sin (dz + a). 



(62) 



(63) 



The phase factor a allows us to consider either a waveguide whose width 
varies sinusoidally 



a = 0, 



(64) 
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or one whose direction changes sinusoidally 

a = t. (65) 

We obtain from equation (54) with 

6 = /3 - ft, (66) 



the Fourier component 



and from equation (55) 



*- = | (67) 



*. = -|..". (68) 

A term of the order a/L <5C 1 was omitted in equations (67) and (68). 
It is apparent that only one spurious mode is excited by the sinusoidal 
wall deflection since condition (66) can be satisfied for only one value 
of m . If condition (66) is not satisfied, <p m and \p m are of the order of 
a/L <SC 1. The fractional power scattered into one spurious guided 
mode due to a sinusoidal wall irregularity is [from equations (48), (53), 
(67) and (68)] 

/AP\ L 2 a 2 k* , 2 t N2 cos 2 K d c o s 2 K m d 2 a ,__.. 

for even modes or [from equations (48) and (60)] 

(*P) U^L {< - i)» cos 2 ^sin 2 ^ gin2 a 



kd + ^ )( M + — 

7o/\ 7m. 

for odd modes. However only one even or one odd mode can be excited 
by one particular sinusoidal wall deviation since it is impossible to 
satisfy the "resonance" condition (66) for more than one mode simul- 
taneously. 

If a = 0, that is if the width of the guide changes sinusoidally, only 
even modes can be excited while sinusoidal deviations from straightness 
(a = ir) couple the even fundamental mode only to odd spurious modes. 
It must also be noticed that for a long period length 

A = | (71) 
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equation (66) can be satisfied only for forward scattering modes. To 
couple to backward scattering modes, the period length D must be 
approximately equal to half the wavelength of the guided modes. 
The fact that only one spurious mode is coupled to the incident mode by 
sinusoidal wall imperfections (it can be shown that the coupling to the 
continuous mode spectrum is also weak if one guided mode can couple 
strongly) allows us to give a much better description of the coupling 
process. 

Since the mode amplitudes C m can change only slowly in the distance 
of one wavelength we can neglect the second derivative of C m in equa- 
tion (34). Labeling the incident mode and the one coupled spurious 
mode 1 we can write the equation system (34) in the following form 

^f = -koxC, , (72) 

' ,tll -*C ll (73) 



dz 



with 



k 2 a , o ,\ cos K Q d cos Kyd /•«! a ,_.> 

«oi = -5" (ttp 1) FT - TTT- „\-u ex P \ l jjj cos %' "§ 



2 



(m + SXm + S. 



The coupling coefficient K 0i of equations (74) holds for coupling from 
an even mode to an even mode 1. The case of coupling from an even 
mode to an odd mode 1 can be treated similarly. In fact, except for 
an unimportant phase factor, we get it from (1/L)[(AP/P) M ]' of equa- 
rion (70). In equation (72) we omitted a term with C on the right-hand 
side, and similarly a term with d was omitted in equation (73). These 
terms would be multiplied by sinusoidally varying functions and would 
describe the local change of phase velocity as the guide dimensions vary. 
These terms give no contribution if we use an average over C and C l 
over the mechanical period length of equation (71). 

Assuming C = \,C\ = at z = the equation system (72) and (73) 
has the solution 

C = cos I koi I z, (75) 

C, = (^Vsin|K 01 |z. (76) 

Total exchange of energy is possible between the two coupled modes. 



3202 THE BELL SYSTEM TECHNICAL JOURNAL, DECEMBER I960 

The distance D over which all the energy is exchanged is given by 

7T 



D = 



& K 01 



(77) 



Finally, we need the power loss to the modes of the continuous spectrum. 
From equations (48), (57), (61), (62), and (63) we obtain 



AP\ a'fc* 2 1V 



cos 2 K t) d 



M + 



To 



i 



2 7 2 « 

COS ad COS - 



+ 



■ 2 . • 2 « 

sin cf/ sin - 



/3(cos 2 crd + h> sin 2 o-rf) /slsin 2 ad + ^ cos 2 a</j_ 
sin 2 [« — (jSo 



*)]£ 



L0 - 0„ - £)] 2 



r/p. 



(78) 



The integration can be performed easily if one realizes that for large 
values of L only a very narrow region in the /3 range near — # — 6 
contributes to the integral. We consider all functions in the integrand 
as constant in this very narrow range and take them out of the integral 
with the exception of 

r 1 2 
Isin [6 - (0 O - 0)] g 



e - (0o - P) 

This remaining integral can easily be performed if we use equation (16) 
to obtain 





dp = -*-dfi. 
p 



Following this procedure yields 



AP 



La 2 k 4 , 2 iv2 cos 2 K n d 



/3 d + 



2 j 2 « 

p COS era COS » 



Tii 



+ 



■ 2 j • 2 a 
psin ao Sin - 



p 2 cos 2 <rd + o- 2 sin 2 o-d p 2 sin 2 ad + a 2 cos 2 ad 



(79) 
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The parameters a and p follow from equations (15) and (16) with 

= O - 0. (80) 

Equation (79) holds only for < fcj we get AP/P = for > fc. The 
most interesting aspect of equation (79) is its linear dependence on L. 
The scattering loss due to the modes of the continuous spectrum acts 
like a true loss process. By contrast, the corresponding equation (69) 
for the loss to guided modes is proportional to L 2 because coupling to a 
guided mode does not result in loss of energy but results in energy ex- 
change between the two coupled modes. Energy loss to one of the guided 
modes is followed by energy gain when the energy exchange has reversed 
itself. 

VI. NUMERICAL EXAMPLES FOR SINUSOIDAL IMPERFECTIONS 

A few numerical examples resulting from equations (74) and (77) 
are listed in Table I. Two different values of the index of refraction n g 
have been assumed, and for each value of the index three different 
values of kd = 2ir(d/\ ) have been chosen so that one, two, or three 
guided modes can exist simultaneously. The mode with O is the lowest 



Table 


I — Numerical Examples 


for Sinusoidal Imperfections 


n„ 


kd 


Pod 


Pxd 


(Srf 


aD 


Remarks 




1.3 


1.729 


— 


— 


— 


Single mode 
operation 


1.5 


1.8 


2.495 


1.916 


— 


6.98 


— 1 coupling 

a = tt 


3.0 


4.336 


3.831 


3.051 
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order even guided mode which is assumed to be incident on the wave- 
guide with sinusoidal wall imperfections. This mode couples to the first 
odd mode with /3j or the next even mode with /3 2 . The values for the 
normalized, dimensionless quantitity (aD)/d 2 [a = amplitude of the 
sinusoidal wall deviation according to equation (62) and (63), d = half 
width of the guide, and D = energy exchange length] have been ob- 
tained with the assumption that equation (66) is satisfied for the two 
modes which are coupled together. Coupling from mode to mode 1 is 
considered only for the case of sinusoidal straightness deviations of the 
waveguide (a = it) while coupling between even modes to 2 is con- 
sidered only for sinusoidal changes of the thickness of the waveguide 
(a = 0). It is immediately apparent from Table I that the energy ex- 
change length D is shorter for a guide with larger values of the refractive 
index. 

To obtain a feeling for the numbers involved in this mode coupling 
phenomenon, let us assume that rc„ = 1.5 and that the free space wave- 
length is Ao = 1m- The value of kd = 1.8 corresponds to d = 0.286^. 
To achieve total exchange of energy between modes and 1 in D = 1 cm 
requires the extremely small amplitude a = 5.72 10~V or a = 0.572 A!* 
The length of the mechanical period in this example is A = 3.1/n. 

Next, let us assume that the index of refraction is n B = 1.01. Using 
again, X = 1/x, we obtain from kd = 15.0 the value d = 2.39m for the 
half width of the waveguide. Requiring again, D = 1 cm, we find a = 
243 A. 

We can look at this problem in a different way. It is unlikely that any 
optical waveguide has a strictly sinusoidal deviation from perfect 
straightness. In fact, the numbers just presented show that it would 
be impossible to produce such a waveguide intentionally. However, we 
have seen [equation (53)] that the mode conversion between two guided 
modes is produced by a Fourier component of the actual deviation func- 
tion. It is therefore not necessary to have a strictly sinusoidal straight- 
ness deviation. Any arbitrary deviation from straightness can be de- 
composed into a Fourier series and the Fourier component at the 
mechanical frequency which satisfies equation (66) is responsible for 
the coupling. In the more general case of arbitrary straightness devia- 
tions, there can be no complete exchange of energy between any two 
modes since power loss to other guided modes and the continuous 



t A mechanical period of a fraction of an Angstrom is somewhat unphysical due 
to the granular nature of matter. However, this result can be restated to say that 
complete power conversion occurs in 0.1 mm if the amplitude is a = 57.2 A. 
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spectrum of modes compete with each other since all of them are coupled 
simultaneously . 

We can now ask the question: What amplitude of the mechanical 
straightness deviation is required to transfer 10 percent of power from 
mode to mode 1 in a distance of L = 1 cm? Again, we use the previous 
examples. From equation (76) [or directly from equations (53) and (77)] 
we obtain 

AP _ | |2 r2 _ 7T^ L 2 _ 

P - I * 01 I u ~ 4 D 2 ' 

For the first example we obtain with n„ = 1.5, AP/P = 0.1, d = 0.286/1, 
and aD/d 2 = 6.98 the value a = 0.115 A.. 1 " This result shows that if 
the Fourier component of the mechanical straightness deviation with a 
period length of 3.1 n is a = 0.12 A (measured over a distance of 1 cm) 
the power loss caused by mode conversion to the first odd mode is 
10 percent. 

For the second example, we use again n„ = 1.01, AP/P = 0.1, d = 
2.39/x, and aD/d 2 = 42.54 and obtain a = 48.8 A. The important 
Fourier component in this case has a period of A = 135/x. The power 
loss to the modes of the continuous spectrum caused by a sinusoidal 
change in thickness of the waveguide (which is very similar to its effect 
as a straightness deviation) can be calculated from equation (79) with 
a = 0. 

Let us consider only one case, n„ = 1.01, kd = 15, A/d = 25. For 
these values we obtain from equation (79) 

4ff -4.(1X10-. 

a L P 

Assuming again AP/P = 0.1 for a guide length L = 1 cm, we obtain 
with d = 2.39/1 

a = 5.46 X 10 -2 m = 546 A. 

This number can be compared to the value a = 48.8 A which gave 
10 percent loss by conversion to one guided mode. However, for a 
meaningful comparison, we must remember that all the Fourier com- 
ponents of a Fourier expansion of the guide imperfections scatter power 
into the modes of the continuous spectrum. The total loss would have 
to be obtained by integrating the scattering loss over the spectral dis- 

* Again it is more reasonable to restate this example to say that 10 percent loss 
occurs over a distance of L = 0.1 mm if a = 12 A. 
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tribution of the Fourier components of the mechanical Fourier spectrum. 
Instead of doing this integration we use a different approach in Section 
VII. 

VII. STATISTICAL TREATMENT OF WALL IMPERFECTIONS 1 " 

Equation (48) gives the relative loss of a guided mode caused by 
a definite (deterministic) distortion of the boundary of a dielectric 
waveguide. A quantity that may be even more interesting is the average 
of equation (48) taken over an ensemble of statistically identical sys- 
tems. 

For simplicity, let us assume that one wall of the waveguide is perfect 
while the other is randomly distorted. If both walls are randomly dis- 
torted, with no correlation between the distortions on opposite walls 
the loss value doubles compared to the case of only one wall being dis- 
torted. If the distortions on opposite sides of the waveguide are per- 
fectly correlated the amount of loss is at most increased four times. 
So to simplify the discussion we assume 

h(z) + d = 0. (81) 

In order to be able to calculate (AP/P) nv , we must evaluate 

(I <p m | 2 >„v = 72 f dz f dz'R{z - I /)«-«''-'-»"' 1 (82) 

Li J Jo 

We assumed that the correlation function 

R(z - z') = <[/(*) - d]W) - d}) av (83) 

depends only on the difference between the coordinates z and z' but not 
on their individual values. 

A change of integration variables allows us to write 

<| <p m | 2 )av = T2 f (L - u)R(u) cos (/3 - i8 m )w du. (84) 

Li J q 

To obtain equation (84) we made use of the fact that R(u) is an even 
function. 

The particular form of R (u) depends on the statistics of the wall 
imperfections. However, all correlation functions have two features in 
common. They all have their maximum value at u = and decrease to 
zero as u — > ■» . If R(u) would not become as u — > °o there would be a 



t An excellent statistical treatment of random coupling effects in metallic wave- 
guides can be found in Ref. 7. 
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systematic distortion of the waveguide boundary instead of the assumed 
random behavior. To get an idea of what one might expect, we assume 
the following form for the correlation function 



(85) 



R(u) = A" «q> I— y 



A is the rms deviation of the wall from perfect straightness and B is 
the correlation length. Using equation (So) we obtain from equation (84) 



;,v L 



OSo - PS + ^ 



1 + 



(00 " ft.)' " £2 



MOSo-fO' + p 



(86) 



where we neglected terms with exp (-L/B) assuming that L/B is 
sufficiently large. In fact if 

L » B, (87) 

equation (86) can be simplified further: 

1 



(I* 



2 . 2£_ 

)nv " BL 



(88) 



(ft, - ft.)' + 



B- 



Using equation (88) we obtain, from equation (53) for the ensemble 
average of the square magnitudes of the even guided modes, 



<|C. 



2\ A. K Jj , 2 i\5 



2B 



cos" K lt d cos K m rf 



(ofc-W + fiXM + S^ + S 



(89) 



The corresponding expression for the odd modes is very similar except 
that cos 2 K m d is replaced by sin 2 K m d and 0„ , K m , and y m are the param- 
eters of the odd modes. 

The total loss caused by coupling to all guided modes supported by 
the dielectric waveguide is the sum over all (\C m | 2 ) av for even as well 
as odd modes traveling in positive (/3 m = +| p„ |) as well as negative 
(0 m = — | p m |) z-direction. It is noteworthy that equation (89) is pro- 
portional to L and not to L 2 . The conversion to spurious guided modes 
by random imperfections appears as a true loss to the incident mode. 

The losses due to the modes of the continuous spectrum are obtained 
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from equations (48), (57), (61), (81) and (88) (with m - 0): 
/AP\ A 2 k*L 2 _ v [" I" p cos 2 K d 

W." «'"--"J-.( a ._ K . + j,)( M+ & 



fcr 



cos 2 <rd , sin 2 «rd 



— 9 ' O 2 ■ 2 "I - » • 2 ) r / tZO « li/Ui 

>,p" cos «rd + cr sin a p sin ad -\- a cos era/. 

The relation between /3, o-, and p is given by equations (15) and (16) 
while 0o , Ko , and t are related by equations (9) and (10) and their 
value is obtained by solving equation (11). The integral in equation (90) 
is extended over /3 from —ktok, the range of real values of the propaga- 
tion constant (in z-direction) of the modes of the continuous spectrum. 
Equation (90) thus includes the losses due to forward as well as back- 
ward scattered radiation. The radiation modes with imaginary values 
of can carry power away from the waveguide only strictly perpendicu- 
lar to its axis. This power loss, if any, is not included in equation (90). 

VIII. NUMERICAL RESULTS FOR THE STATISTICAL CASE 

Figures 4 through 9 show numerical evaluations of equations (89) 
and (90). These figures can be grouped into two classes. Figures 4 
through 6 are drawn for a dielectric waveguide whose index of refraction 
is n a = 1.01. Figures 7 through 9 apply to a waveguide with n„ = 1.5. 
Within each of these two classes, the kd value was chosen to allow for 
three different cases. Figures 4 and 7 apply to waveguides which can 
support only the lowest order guided mode. In this case there is power 
lost only to the modes of the continuous spectrum. Figures 5 and 8 
apply to waveguides supporting two guided modes and Figs. 6 and 9 
apply to waveguides supporting three guided modes. Each figure shows 
the normalized loss caused by scattering into modes of the continuous 
spectrum as solid lines and the loss to the possible guided modes as 
dotted lines. Also shown are the ratios of backward to forward scattered 
power as solid lines for the modes of the continuum and as dotted lined 
for the guided modes. The total power lost to the lowest order guided 
modes is the sum of the losses to the continuum and the spurious guided 
modes. 

Several remarkable features of these loss curves are worthy of a com- 
ment. The losses caused by the modes of the continuum as well as by 
the guided modes peak at certain values of the correlation length B. 
The location of these peaks are different, however, for the continuum 
and guided modes. 
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Fig. 4 — Normalized radiation loss (d*/A 2 L) (AP/P) and ratio of backward to 
forward scattered power AP'/AP* as functions of the normalized correlation 
length B/d for n„ = 1.01 and kd = 8.0. Single guided mode operation (d = half 
width of waveguide, A = rms deviation of one waveguide wall, L = Length of 
waveguide section, n g = index of refraction of waveguide, k = free space 
propagation constant). 

The losses to the guided modes increase with increasing number of 
guided modes supported by the waveguide. However, the losses caused 
by the continuum of modes also increase as an increasing number of 
guided modes can be supported. This increase is less rapid, however, 
as one might expect because of the dependence of equation (90) on the 
fourth power of k. The fourth power dependence on frequency (or inverse 
wavelength) is typical for Rayleigh scattering by small particles, and 
it is not surprising that we encounter it here. 

Finally, it is apparent from the curves showing the ratio of back- 
scattered to forward scattered power that forward scattering is pre- 
dominant for large values of the correlation length. The ratio of 
AP~/AP + levels off for large values of B. In some of the curves the 
leveling of the AP~/AP + curves occurs out of the diagram but it is a 
common feature of all the curves. For small values of the correlation 
length there is as much scattering in the forward as in the backward 
direction. 

For many practical applications, a waveguide supporting only one 
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Fig. 5 — Normalized power loss and ratio of foreward to backward scattered 
power for radiation (solid curves) and spurious guided modes (dashed curves). 
Two guided modes (n„ = 1.01, kd = 15). 
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Fig. 6 — Similar to Fig. 5. Three guided modes (n, = 1.01, kd = 23). 
two guided mode loss; continuum loss. 
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Fig. 7 — Similar to Fig. 4. One guided mode (n, = 15, kd = 1.3). 
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Fig. 8 — Similar to Fig. 5. Two guided modes (n„ = 1.5, kd = 1.8). 
one guided mode loss; continuum loss. 
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Fig. 9 — Similar to Fig. 5. Three guided modes (n g = 1.5, kd = 3). 
two guided mode loss ; continuum loss. 

guided mode may be of most interest. Let us assume X = 1m- Figure 4, 
holding for kd = 8.0 and n g = 1.01, applies to a waveguide whose half 
width is d = 1.27 n. Taking the worst possible case of B/cl = 9 or B = 
11.4/x, we find from Fig. 4 

d 3 AP 3 

If we want to know how much rms deviation A of one wall of the guide 
would be required to cause a 10 percent loss (AP/P = 0.1) in one centi- 
meter of waveguide (L = 1 cm) we find A = 5.85 X 10~ 2 (t = 585 A. 
The ratio of A over d gives an idea of the relative tolerance require- 
ments: 



j = 4.6 X 10 -2 = 4.6%. 

If the waveguide were to conform to the conditions of Fig. 6, we would 
have for X = 1/x a half width d = 3.66/z. The losses caused by the two 
spurious modes are of the same order of magnitude as the radiation losses 
caused by the continuous spectrum. For B/d = 10 or B = 36.6/x we get 
a total loss of 
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d * £ = 3.4 X 10". 



A*L P 
To cause AP/P = 0.1 for L = 1 cm requires that 

A = 1.2 X 10" V or ^y = 3.28%. 

The relative tolerance requirements are, therefore, approximately the 
same in both examples. 

As a last example let us use Fig. 9 corresponding (X = 1m) to a wave- 
guide with n = 1.5 and a half width d = 0.477^. For B/d = 1.3 or 
B = 6.2/i we find for the total loss 

d ' AP = 2.3 X 10-. 



A'L P 
We get AP/P = 0.1 with L = 1 cm for 

A - 2.18 X 10~V = 21.8 A or ^ = 0.457%. 

The perturbation theory, strictly speaking, holds only for small values 
of AP/P. However, it is reasonable to expect that the power scattered 
into the radiation modes escapes sufficiently rapidly so that no appreci- 
able amount of power reconversion from the radiation field to the guided 
mode occurs. The incremental power loss, AP/P = — aL, is therefore 
the same for any section of the guide so that we obtain the total scatter- 
ing loss into the continuum of radiation modes P = P l) e~ aL . We may 
now ask how much rms deviation is required to cause a radiation loss 
of 10 dB/km or« = 2.3 km -1 = 2.3 X 10~ 5 cm" 1 . Using B/d = 10, 
corresponding to the top of the loss curve of Fig. 4, we obtain the 
equation 

42 X 2.3 X 10" 5 = 6X 10" 3 
A 

so that (X = 1/x, n„ = 1.01, kd = 8.0, d = 1.27 X 10~ 4 cm) 

4 = 6.98 X 10" 4 or A = 8.86 X 10" 8 cm = 8.86 A. 
d 

This figure dramatizes the stringent tolerance requirements of dielectric 
waveguides for long distance optical communications. In fact, such 
tolerances seem impossible to obtain. One can only hope that the cor- 
relation length can be kept far from the worst possible value of B/d = 10 
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(in this example) so that these extremely stringent tolerance require- 
ments might be eased. 

IX. CONCLUSION 

We have analyzed the losses suffered by the lowest order symmetric 
mode propagating on a dielectric slab waveguide caused by imperfec- 
tions of the waveguide boundaries. The analysis was simplified by assum- 
ing that there is no change in either the dielectric slab or the guided and 
unguided fields in one direction parallel to the slab. This assumption 
causes all our conclusions to be optimistic since variation of the slab 
in this direction can only cause additional losses. However, we expect 
that the results of this analysis give at least the correct order of magni- 
tude of the actual scattering losses. 

The statistical analysis was limited to a study of the effects which 
an exponential correlation function might have on the waveguide losses. 
The actual form of the correlation function may be quite different from 
this assumed exponential shape. 1 " Conclusions regarding loss predictions 
are further hampered by a lack of knowledge of the expected correlation 
length. 

However, our analysis does lead one to conclude that scattering losses 
suffered by optical fibers or other dielectric waveguide structures may 
be very serious. Deviations of the waveguide wall in the order of a few 
percent can cause a power loss of 10 percent or 0.46 dB/cm if the wall 
imperfection can be described by an exponential correlation function 
with a correlation length to guide half width ratio of approximately 
B/d = 10. An rms deviation of A = 9 A causes a radiation loss of 10 
dB/km if the free space wavelength is X = 1m and the guide has an 
index of refraction of n B = 1.01 (with vacuum on the outside). The 
width of the slab in this last example is 2d = 2.54^. 

The mode coupling and radiation loss theory has been experimentally 
confirmed at microwave frequencies. A report on these measurements 
is given in Ref. 8. 



t Several other correlation functions have been tried and it was found that the 
results are insensitive to the particular choice of the function for values of B/d 
less than the value corresponding to the loss peak. In particular, the maximum 
loss value and the position of this loss peak were the same for different correla- 
tion functions. However, the loss values for B/d larger than the value correspond- 
ing to the maximum of the curve are very strongly dependent on the choice of 
the correlation function. 
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